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1 Introduction 

A great number of disordered systems exhibit very long relaxation times as 
some critical temperature is approached. Below this temperature equilibrium 
is no longer reached within finite time and the behaviour of such systems 
becomes non ergodic. Among those are diluted alloys of magnetic ions in a 
non magnetic matrix, so called spin-glasses P, supercooled liquids |2] enter- 
ing some glassy state at low temperature or particles moving in a random 
potential exhibiting a transition from diffusion to creep or pinning [3 E] . The 
temporal behaviour of such systems is often referred to as glassy dynamics. 
Below the characteristic temperature aging is observed If, for instance, a 
spin-glass is cooled in a magnetic field and kept for some waiting time, the 
complete decay of the magnetisation after the field has been switched off, is 
hindered for times of the order of the waiting time. Similar phenomena are 
observed for the deformation of glasses under the influence of applied forces. 

Glassy dynamics can also be of relevance for systems and phenomena 
outside physics. For some combinatorial optimisation problems simulated 
annealing is used in order to flnd good solutions. This means that some 
stochastic dynamics, characterized by some "temperature", is used in order 
to flnd minima of a properly deflned cost function. If such a system exhibits 
a freezing transition, it is not effective to spend much time in simulated 
annealing below this temperature and investigating the dynamics slightly 
above the freezing transition is much more efficient Other examples of 
systems of interest outside physics are models of neural networks or certain 
aspects of markets and other economical systems [7j. 

The systems mentioned above are classical. There is also interest in quan- 
tum mechanical disordered systems, for instance in certain spin-glasses or 
in interacting tunnelling systems being of relevance for glasses at very low 
temperatures. This lecture deals, however, with classical dynamical systems 
only. 

A breakthrough of our understanding of the low temperature properties of 
glassy systems came from replica theory and the replica symmetry breaking 
scheme proposed by Parisi [8^. This theory focuses on the evaluation of the 
free energy averaged over some frozen disorder. The picture emerging from 
this theory is a decomposition of phase space into so called pure states or 
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ergodic components, separated by barriers of infinite hight in the thermo- 
dynamic hmit. Investigating the overlap among those states, an ultrametric 
organization reveals. 

An alternative treatment of systems with quenched disorder is via dyn- 
amics jni Uni im II2I ■ THs win be the subject of this lecture. As it turns out 
there are many features and results common to both treatments, but there 
are also differences which will be discussed. 

In the first part I present some of the systems of interest, in particular spin- 
glasses, supercooled liquids and glasses, drift, creep and pinning of a particle 
in a random potential, neural networks, graph partitioning as an example 
of combinatorial optimisation, the K-sat problem, a prototype problem from 
computer science, and finally the minority game as a model for the behaviour 
of agents on markets. 

The second part deals with the dynamics of the spherical p-spin-glass with 
long ranged interactions. This model is a prototype for glassy dynamics. 
The equations of motion for correlation- and response-functions are derived 
and equilibrium solutions above the freezing temperature are investigated. 
Below the freezing temperature the system behaves for short time as if it were 
in equilibrium, for longer time its off equilibrium nature becomes apparent 
and correlation- and response-functions depend on waiting time. This is the 
manifestation of aging. 

At the end of this part I discuss the so called crossover region and aging in 
glasses. Within mode coupling theory, which is discussed in the contribution 
of R. Schilling in this volume, and which shows great similarity to the 
p-spin model, a dynamical transition is found at some critical temperature Tc. 
This is sometimes referred to as "ideal glass transition". Comparison with 
experiments shows that this transition is actually smeared out and the fitted 
Tc is well above the actual glass transition region. Some smearing out also 
results from an extended version of the mode coupling theory, which applies, 
however, to equilibrium only. I am going to discus a modified version of the 
p-spin model, which is able to take account for this rounding as well as aging. 



2 Examples of Disordered Systems 
2.1 Spin-glasses 

Typical spin-glasses are diluted alloys of magnetic ions in a non magnetic 
matrix, e.g. Agi^xMux with x ~ 3%. The magnetic ions interact via the 
RKKY-interaction which is mediated through the conduction electrons. It is 
of the form 



shown in FigO For random distances between pairs of magnetic ions, 
their exchange interactions Jij = J{rij) are also random variables including 
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positive and negative values. This leads to frustration, as shown in the insert 
of FigHlfor a triple ijk of ions with Jjj > 0, Ji^k > and Jj^k < 0. 




Figure 1: RKKY-interaction 

At high temperatures spin-glasses typically show a Curie-law for the mag- 
netic susceptibility x ~ '^/T, as shown in FigI21 for CuMux. Below some 
freezing temperature Tg the system is no longer in equilibrium. If the system 
is cooled in a small field, the magnetisation M = XfcB stays more or less 
constant (a and c in Fig|21). The ratio xfc = M/B is called field cooled 
susceptibility. 
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Figure 2: Susceptibility of CuMn^ for different values of the concentration 
X of the magnetic ions. 

If, however, the system is cooled without applied field, and the field is ap- 
plied only after some waiting time t^, the resulting magnetisation is reduced, 
i.e. the zero field cooled susceptibility Xzfc = M/B < xfc ip and d in 
Fig|21). Approaching the critical temperature Tg critical slowing down can be 
observed, see FiglHl 

The off equilibrium character of the low temperature phase is most clearly 
demonstrated by the decay of the remnant magnetisation. In this experiment 
the sample is rapidly cooled in a magnetic field to a temperature T < Tg. 
After some waiting time tw the field is switched off, and the decay of the 
magnetisation is observed as function of time t elapsed after removal of the 
field. Following some rapid initial decay, not shown in FigU the magnetisa- 
tion stays almost constant up to time t ~ t^. 
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Figure 3: Real part of the susceptibility x'{^) ^oi Feo.sMNo.sTiOs ^ISj. The 
critical temperature is Tg = 20.7 [K] 
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Figure 4: Decay of the remnant magnetisation for various waiting times 
indicated by arrows [TB] . 

Investigating more elaborate temperature programs T(t) a variety of other 
striking memory and aging effects are observed jT7j. Similar aging phenomena 
are also found in glasses 

Theoretical investigations are typically based on Ising models with random 
interactions. Their energy is 

= -i ^ JijCTiaj ai = ±1. (2) 

The interactions Jij are assumed to be gaussian distributed random variables 
with 

Ji,j = Ji,jJk,l = \{Si,kSj,l + Si,lSj,k}Wij. (3) 

Calculating e.g. the free energy 




F = -kBT In J2 e-^^) (4) 

{a=±l} 

the problem arises, that the logarithm of the partition function has to be 
averaged over the gaussian disorder. Similar problems exist in evaluating 
expectation values, e.g. 

This can be done using the replica trick 3J. 

As an alternative, one can examine dynamics in the form of stochastic 
processes having the Boltzmann distribution ~ e~^^ as stationary solution. 
For an Ising model, Glauber dynamics can be used. It is is given by a single 
spin-flip master-equation such that 

A = (^tanh (/?^ J,,a,(t))^ . (6) 

Using a path integral representation of this process, the average over the 
stochastic interactions can easily be performed. This will be subject of the 
second part of this lecture. 



2.2 Supercooled Liquids and Glasses 

If a liquid is cooled sufficiently slow, it may avoid crystallisation and enter 
the state of a supercooled hquid. With decreasing temperature the shear 
viscosity r]{T) increases and reaches a value of i]{Tg) ^ 10^^ [Pa sec] defining 
the glass temperature Tg. At this value, plastic fiow can hardly be observed in 
laboratory experiments. In some glasses, so called strong glasses, the viscosity 
follows more or less an Arrhenius law //(T) ~ e^^°- ^Hl- A typical strong glass 
formers is Si02 with an activation energy Ea ~ 4eV corresponding to the 
binding energy of a covalent bond. Fragile glass formers, on the other hand, 
show pronounced deviations from this law. Typical examples are organic 
molecular glasses, ionic glasses, polymers or proteins. Some examples are 
shown in Fig|3 

There are other attempts to define a glass transition temperature, e.g. 
identifying a rapid drop in the specific heat or thermal expansion [20] • Ac- 
tually the glass temperature defined in one way or another usually depends 
on the cooling rate, and it is not clear at all, if some finite Tg exists in real 
glasses. There is no need to discuss this further in this lecture, since theo- 
ries of the structural glass transitions are covered in the contribution of R. 
Schilling ^3] in this volume. 

There is, however, one aspect of interest in the present context. This is the 
"Ideal Glass Transition" found in mode coupling theory (MCT) ^31^. The 
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Figure 5: Angell plot: Logarithm of the viscosity versus inverse temperature 
|18j . A straight hue corresponds to Arrhenius behaviour. The fat symbols 
indicate the MCT-crossover temperature |19j . 



resulting equations resemble those found for the dynamics of spherical p-spin 
interaction spin-glass ^2^2 be discussed later in this lecture. This theory 
yields diverging time scales as some critical temperature Tc is approached, and 
non ergodic behaviour for T < Tc. 

An onset of critical slowing down, usually referred to as a-relaxation, can 
be observed in fragile glasses, e.g. in the dielectric susceptibility x"{^) of the 
ionic glass CKN ||^ shown in FiglHla. Fitting the critical temperature 
and other parameters of mode coupling theory, as indicated in FiglHlb, yields 
the critical temperatures shown in Fig|21 
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Figure 6: (a) Imaginary part of the dielectric susceptibility yl^'iiS) of the ionic 
glass CKN j2Sl- The a-relaxation shows up as the broad peak at the lower 
frequency, (b) The frequencies of the maximum cUmax and minimum tUmin of 
fitted to cUmax ~ (T - TJ'' and uj^,^ ~ (T - T^)'''. 



Actually an extended version of the mode coupling theory 121] yields 
a rounding of the singularity near T^. and equilibrium behaviour for temper- 
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atures below Tc, in contrast to the findings of tlie simplified mode coupling 
theory or the p-spin-glass. On the other hand, aging can be observed in 
glasses as well, indicating off equilibrium properties at low temperatures. An 
example is shown in FigjTI We shall come back to this point later. 
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Figure 7: Aging observed in PVC A shear stress is applied after some 
waiting time following a quench from the liquid state. Then the resulting 
deformation is observed as function of time t. 

Assuming pair interactions, the (potential) energy may be written as 

H = ^ J drdr' V{f— r') n{r) n{r') (7) 

where n(r) is the density at point r. 

If one is interested in slow dynamics only, it is sufficient to investigate the 
limit of overdamped motion, i.e. the Langevin equation 

—n{f, t) = Vn(r, ■ V / dr' V{r — r') n{r', t) + fj{r, t) ■ Vn{f, t) (8) 

with fluctuating forces such that 

(r/(f,t)) = {f]{f,t) -f^ifX)) = '^T5{r -r)5{t-t'). (9) 

Investigating time dependent correlation-functions, the nonlinear term in the 
above equation is treated via mode coupling theory. As mentioned, the result- 
ing equations of motions are identical to those obtained for the p-spin-glass 
(in equilibrium). This is remarkable, because the Hamiltonian of the spin- 
glass has built in quenched disorder, whereas the Hamiltonian of the glass 
does not, and the disorder appears to be self generated. 

The relevance of the findings of mode coupling for a glass transition at Tg 
is not completely obvious. Other concepts have been proposed to understand 
the glass transition, e.g. activated transitions among "inherent states" |2HI 
or models with kinetically constraint dynamics j2Zl I2H] ■ 
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2.3 Motion of a particle (manifold) in a random poten- 
tial 



Another example for glassy dynamics is the motion of a particle in a correlated 
random potential [31I3|. Using again the limit of strong damping, a Langevin 
equation may be used 

v{t) = ^f(t) = -VV{f{t)) + K + ff{t) (10) 
where V{r) is assumed to be a gaussian random potential with moments 



V{r) = nrlV{P) = j^^-^^, (11) 

— * 

and K is some external force pulling the particle. The thermal noise is again 
given by Eq.Q. 




Figure 8: Drift, creep and pinning for a particle moving in a random potential. 

Depending on temperature T and the exponent 7, characterising the range 
of the random potential, various types of motion are found Q , see FiglSl For 
T > 7^(7) the average velocity v{K) ~ K, i.e. there is a finite friction 
constant. For T < Tc{'y) and 7 < 7c the particle moves only if the force 
exceeds some critical value K > -ft'pinning(7', 7). For T < Tc(7) and 7 > 7c 
creep is found, i.e. v{K) ~ with 77 < 1. 

A similar freezing transition shows up for a chain or some other manifold 
moving in a random potential |29j . 

2.4 Neural Networks 

Considerable effort has been put into the study of formal neural networks 
dealing with random patterns e.g. the Hopfield model jSI] as a prototype 
of an associative memory. 

Representing the "on" and "off" state of a neuron by the two states of an 
Ising spin a, = ±1, the energy is assumed to be that of an Ising spin-glass 
Eq.Q. The couplings Jij are, however, determined by the patterns to be 
memorised. Let = ±1 be a set of random patterns, where i = 1- --N 
labels the neurons and fi = 1 ■ ■ ■ A the patterns. Then the choice of the 
couplings Jij such that 
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ensures, that a fixed point of tlie dynamics, defined by Eq.®, exists near 
each pattern /x at least as long as A < with ~ 0.13. Each fixed point 
is surrounded by some basin of attraction, i.e. if the initial state contains a 
sufficient part of a pattern, the complete pattern is reconstructed with a small 
number of errors. The size of the basin of attraction depends on the loading 
a = A/N [32]. For a > ac the memory is overloaded and looses its function 
completely. The dynamics becomes glassy, which is not too surprising, since 
the couplings Jij can be positive as well as negative, resulting in frustration. 

Methods developed for spin-glasses, replica theory and dynamics, have 
also been used to investigate this and other types of neural networks. In 
particular bounds of learning have been determined. 

2.5 Combinatorial Optimisation Problems 

Combinatorial optimisation problems show up in various technical applica- 
tions. The problem is, to find the "best" out of a discrete set of states. 
Finding the ground state of an Ising spin-glass or solving the travelling sales 
man problem, i.e. finding the shortest path connecting a set of towns, are ex- 
amples. In many cases the problem is A^P-complete, i.e. the optimal solution 
can not be found with computational effort growing with the size N of the 
problem to some power. This means that for larger systems it is virtually im- 
possible to find the optimal solution. In many applications it might, however, 
be sufficient to find good solutions in polynomial time. One of the standard 
procedures is simulated annealing It is implemented by constructing a 
cost function (energy), and using some stochastic dynamics, characterised by 
a temperature. This temperature is slowly decreased, and hopefully states 
with low values of the cost function are found. 

Actually such a system might undergo a transition to glassy dynamics 
below some critical temperature Tc. If this is the case, extending simulated 
annealing to T < Tc is not efficient, and searching for good solutions at T ^ Tc 
is more effective 6j. 

A standard problem in this context is graph bipartitioning [34j. Assume 
electronic components i = 1 ■ ■ ■ N have to be placed on two chips A and B. 
The aim is to place them such that the number of connections between A and 
B is minimal. A cost function for this problem can be defined by mapping it 
to an Ising system with 



cTj = 1 if i on A 



ai = —1 if i on B. 



(13) 



and interactions 



Jij = 1 if 2 and j are connected 



Jj j = else. 



(14) 



The resulting cost function is 




(15) 
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For the process of simulated annealing again Glauber dynamics, Eq.®, may 
be used. 

2.6 Random K-sat Problem 

Another standard A^P-problem is the i^'-sat problem [33] . This problem deals 
with N Boolean variables represented by cXi = ±1. There are A clauses with 
K literals per clause among the Boolean variables or their negations, e.g. for 
K = 3: R = {(Ti OR (jj OR NOT cr^ }. A clause R^{{a}) can be represented 
as 

{1 if cTj in Rf, 
-1 if NOT a, in R^ (16) 
else . 

N 

The clause R^{{a}) is not fulfilled if = K, and a cost function 

i=l 

AN ^2 

/i=i 1=1 /I 1=1 

may be used. Finding solutions of the i^-sat problem, again simulated an- 
nealing can be used, and finding a ground state with H = means that the 
corresponding choice of the Boolean variables satisfies all clauses. The formal 
similarity to the Hopfield model is obvious. For random clauses a critical 
ac{K) exists, such that for a = A/K < ac{K) solutions can typically be 
found by the polynomial simulated annealing procedure. 

2.7 Minority game 

Methods, developed for disordered systems, have also been applied to various 
problems in economy, especially concerning financial markets. Activities of 
this kind are usually subsummized under the notion econophysics. As an ex- 
ample the minority game j7j is discussed. This is a repeated game, where N 
agents have to decide which of two actions, e.g. buy or sell, to take. The goal 
is to be in the minority. Each agent knows the outcome of the last L games 
and can select among two strategies, where the set of strategies is different 
for each agent. For each new game the agent selects the strategy which would 
have been most successful over the last L games. Again this problem can be 
mapped onto the dynamics of an Ising spin system with quenched disorder, 
the disorder being the random set of strategies for each agent. Depending on 
N and L freezing transitions are found. 
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3 Dynamics of the ^^-Spin Interaction Spin- 
Glass 



3.1 Soft Spin Models 

In the previous part of this lecture we have seen many different examples 
of disordered systems exhibiting transitions to glassy dynamics and off equi- 
librium behaviour. Some of these systems are based on Glauber dynamics 
with transition probabilities among a discrete set of states. Glauber dynam- 
ics is very appropriate for Monte-Carlo simulations, but difficult for analytic 
investigations In the following we therefore investigate a model based 
on continuous degrees of freedom. In equilibrium the resulting equations of 
motion are identical to those derived from mode coupling theory for glasses 
O Cni • The model is, however, more general in the sense that it allows to 
study off equilibrium properties and aging 0. 

Actually we investigate a whole class of models. The degrees of freedom 
are scalar real variables yjj representing soft spins or other variables. The 
energy (Hamiltonian) is 

H = ^U {ipi) - ^ hiifi Jij^i^j - ![ X] Jijk^i^j^k (18) 

i i i,j i,j,k 

where U{ip) is some local potential constraining the distribution of the (pi. 
Ising spins for instance can be approximated by choosing U{ip) as deep double 
well potential with minima at = ibl, e.g. U{ip) = U [1 — </5^)^. For a 
spherical model U{(p) = |/iV^^ with /i such that {ip"^) = 1. 

We allow for multi spin interactions Jijk ■ ■ ■ involving two and more spins 
or particles. The interactions are Gaussian random variables with 

Jij = JijJkl = 2^^hk^j,l + ^i,l^j,k}W^j 

Jijk = JijkJlmn = ^{Si,lSj^rnSk,n + ■ ■ ■} (19) 

For real spin- glasses the range of the interactions is finite, and one may choose 
W}p = W if site i and j are next neighbours and w!:^^ = else (Edwards- 
Anderson model [SE])- For the SK-model (Sherrington-Kirkpatrick [HZj) on 
the other hand interactions with unlimited range are used. Models with long 
ranged interactions (or models in infinite dimensions) have the advantage 
that mean field theory yields the exact solution. For other problems like 
neural networks or combinatorial optimisation problems space has no meaning 
and interactions among all elements are appropriate. The strength of the 
interactions has to scale with the number A^ of elements in order to allow for 
a non trivial thermodynamic limit N ^ oo 

W^^ = N~P/^ Wp. (20) 

Models of this kind are often referred to as mean field models. 
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3.2 Replica Theory 

Equilibrium properties can be derived from the free energy averaged over the 
disorder 

F{T, h) = -kBT\n Z{T, h, J)\ (21) 

The problem is to perform the average of the logarithm of the partition func- 
tion instead of averaging the partition function itself and taking the logarithm 
afterwards. The replica trick jHl ISZ| uses the identity 

ln(a;) = lim fa;" - iV (22) 

For integer n the partition function may be computed by replicating the 
system n times 

n „ 

= JJ / d<---dv7^e-^^"^(^"). (23) 

a=l ^ 

This quantity can now be averaged easily over the gaussian disorder. Since the 
disorder is the same for each of the replica, the resulting effective Hamiltonian 



a a, (3 

couples different replica. The partition function can be evaluated for mean 
field models introducing the order parameters 

q^,p = (^^f,^^)\ (25) 

The problem is symmetric with respect to permutations of the replica. At 
high temperatures T > the order parameters are symmetric with respect 
to permutations as well, i.e. 

qa,/3 = qoSa,f3 + ? (^1 " (^a./s) • (26) 

At low temperatures T < two different replica symmetry breaking schemes 
have been proposed |H], a one step symmetry breaking and alternatively a 
hierarchical symmetry breaking in p steps with p — > oo, the so called full 
symmetry breaking. At the end the limit n — has to be taken extrapolat- 
ing the results obtained originally for integer n to real n. The results for the 
low temperature state with full symmetry breaking are interpreted in terms 
of a state composed of distinct pure states with ultrametric organisation with 
respect to their overlap. Pure states are interpreted as regions of phase space, 
valleys, surrounded by barriers of infinite hight. Within each valley the sys- 
tem is in equilibrium. A corresponding interpretation is obtained from the 
following investigation of dynamics. For further details on the replica theory 
and results see jHl CHI ■ 
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3.3 Langevin Dynamics and Path Integrals 

An investigation of dynamics yields additional information about the tempo- 
ral behaviour of correlations and especially about the critical slowing down at 
the freezing transition. Below the transition temperature T < Tc signatures 
of non equilibrium are expected, in particular aging. 

Specifying a Hamiltonian for a classical system does not specify its dy- 
namics uniquely. The dynamics has, however, to be chosen such that the 
equilibrium state ~ e~^^ is stationary. The simplest form is a Langevin 
equation corresponding to the overdamped motion of particles or spins 

+ (27, 

where the friction constant has been put to 1. The fluctuating forces, repre- 
senting a heat bath, are assumed to be Gaussian distributed random variables 
with 

{r]i{t)) = {v^{t)vAt')) ='^T6^^j 6{t - t'). (28) 

In contrast to the quenched disorder of the interactions J..., the fluctuating 
forces r]i{t) depend on time. Averages over the J... and the thermal motion, 
are distinguished by different notations, ^^'^ and (■■■), respectively. 

It is convenient, especially in view of the average over the quenched dis- 
order, to introduce a path integral representation [21] • In the following the 
derivation is outlined, dropping the site index i for a moment. Let 

^{t)=^it;{r]y,^o) (29) 

be solution of Eq.(j2Zj) with initial condition ip{to) = (fo- The thermal average 
of a product of (^-variables at different time is then 

{vitMt') ■■■)=j d^Mvo) jjv{ri] e^^'^^^^^Vlt; {^}; v^o) ^{t': {^}; y^o) ■ • • 

(30) 

with 

>V(W;^o) = -4^ / dtriitf. (31) 

Po(v'o) is some distribution of initial values. 

Instead of integrating over the fluctuating forces Ti{t), it would be more 
convenient to have a path integral over the dynamical variables ^p{t). This 
can be done by introducing imaginary auxiliary variables <^(t), and then 



{^{mt') ■■■)=] ^VoP.ivo) j ^, V} e^«^'^'''>^*° V W m ■ ■ ■ (32) 

with 

W{{^, ^, r/}; to) = - ^ dt { ^r^itf + m [m - F(^(t)) - r/(t)] } . (33) 
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The path integral over the auxihary imaginary variables 0{t) can be viewed 
as a product of (5-functions insuring the fulfilment of the equation of motion, 
Eq-dZZD, at all time. 

The final step is to integrate over the fluctuating forces r]{t), which can 
easily be done by completing the square in Eq. ()33|) . and now 



to 



{^{tMt') ■■■)=/ d^oPo(^o) / / e^«^'^>^*o)^(t) m ■ ■ ■ (34) 

with 



W({^, to) = dt\ T^{tf - m m - Fi^it)) . (35) 



to 

The path integral extends over real functions (p{t) with initial condition 
(p{to) = (po, and imaginary functions ip(t) with unrestricted initial condi- 
tions. The path integral goes over functions in the interval to < t < ti with 
ti greater than the latest time argument in the expectation value Eq. ()34|) . 

The auxiliary fields have a physical interpretation. Allowing for a time 
dependent external field h(t) in Eq. (fTH|) . response functions can be expressed 
as expectation values of the original and auxiliary fields, e.g. 

(36) 

The expectation value on the right hand side is of the form Eq. ()34|) . The 
auxiliary fields are therefore denoted as response fields. 

Response functions have to be causal, i.e. Eq. ()36|l has to vanish for t' > t. 
This is in accordance with the Ito calculus for the stochastic equation of 
motion Eq. ljTfjl . assuming that the action of the forces is retarded [SHI 122] • 

The formulation given above is completely general in the sense that it is 
not restricted to equilibrium or small deviations from equilibrium. In general 
correlation- and response-functions are not related. In equilibrium, however, 
they have to obey fiuctuation-dissipation-theorems (FDT) 

^±q{t,t')=T{tX) (37) 

where correlation- and response-function are defined as 

q{t, t') = {ip{t)ip{t')) r{t, f) = ivimt')) . (38) 

A sketch of the proof is as follows: The first step is to show that for 
equilibrium initial condition 

Poiifo) = Z-ie-^^(^«) (39) 

the initial time is arbitrary as long as it is earlier than the earliest time in 
expectation values of the type Eq. ljH^ or Eq. (jHU|) . Assume that the initial 
distribution at some time < to is given by Eq.p^. For tQ < t < to one 
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replaces 0{t) ^ ip(t) + l3ip{t) in the path integral Eq. ljH^ and in the action 
Eq. (jH3|) . This results in 

W({^,^};to)->V({^,^};g = rdt{Tm^ + m\m+nvit)) 



- /3Hicp{to))+f3Hicp{t',)). (40) 

The last two terms in Eq. ()40|) originate from integrating a contribution 
j3(f(t)F{ip(t)) = —[3dH{ip{t))/dt which appears if the above substitution 
is performed. Evaluating the path integral over ip{t) or actually ip{t) for 
I'q < t < the only contributions come from (p{t) = 0. The last two terms 
in Eq. ()40|) replace the initial equilibrium condition Eq. ()39p at tj, by the cor- 
responding one at to. If the system is ergodic, the limit ~^ t^^Y be 
taken and, as a consequence, the initial condition becomes irrelevant. 

Consider now a small change of the external field 6h{t) = 6hQ for t < to 
and Sh{t) = for t > tg- Assuming equihbrium 

rdt'r{t,t')=Pq{t,to). (41) 

For the first expression the initial time t'^ — —oo is used. The second is 
obtained by differentiating the equilibrium initial condition at to with respect 
to h. The FDT is obtained by differentiation both expressions of Eq. (j4H) with 
respect to to. 

Deriving this result no explicit time dependence of the Hamiltonian is 
allowed and the system has to be ergodic. 

3.4 Average over Quenched Disorder 

As pointed out already in the context of the replica theory, the calculation 
simplifies considerably for mean field models with long ranged interactions, 
specified in Eq.fl2U|). In the following we investigate models of this kind 

IHini Ill- 
Provided the initial condition does not depend on the disorder, the average 
over the interactions J^^^ can easily be performed using the general relation 



valid for Gaussian distributions. This yields, for the Hamiltonian Eq. ()18|1 and 
the disorder specified in Eq. lfT^ and Eq. (p]j) . the effective action 

w({^,^})' = [dtj2{THt?-Ht)[m-h-u'{v,{t))]] 

+ j^tj\t'Y,{mw'[j^Y.^,{t)^,{t'))ut') (43) 
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with 

w{x) = J2¥ ^^^^ (^^) 

p 

Correlation- and response- functions calculated with this action are now dis- 
order averaged quantities 



3.5 Dynamic Mean Field Theory for Spherical Models 

Computing local correlation- and response-function, i.e. averages involving 0i 
and (pi at a single site i only, a saddle point evaluation of the corresponding 
path integral is possible. The local time dependent functions are obtained 
from an effective single site action 



Weff({(^i,(/?i}) = Jdt\T^iity -ipiit)^ipiit)-h-U'{ipiit)) 

+ Jdtj'dt'[0,{t)W'{q{t,t')') 0,{t') (46) 



+0,{t)r{t,t')W"(q{t,t'))^^{t' 



which has to be determined selfconsistently. This is exact for mean field 
models in the thermodynamic limit N —>■ oo. 

Compared to conventional mean field theory, e.g. for magnets, the present 
theory is much richer because the order parameters are the correlation- and 
response-functions and therefore functions of two time arguments. 

For spherical models with U{x) = the effective action is quadratic, 
which simplifies the calculation further. In particular it allows to write down 
the resulting dynamical mean field equations in the following closed form 



d 

(— + /i(t)) g(t,t') = hm{t') + dsK{t,s)q{s,t') 

+ ds[M{t,s)r{t',s) + K{t,s)q{t',s)^ (47) 

A ft 

{^^+Kt))rit,t') = j^^dsK{Us)ris,t') (48) 

+ /i(t)) m{t) = h + dsK{t,s)m{s) (49) 



with 



K{t,t') = W"{q{t,t'))rit,t') 

M{t,t') = W'{q{t,t')). (50) 
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The spherical constraint q{t,t) = 1 yields 



fi{t) = hm{t) + T + ds s)q{t, s) + M(t, s)r (t, s) | . (51) 

The above dynamical mean field equations are a set of coupled non lin- 
ear integro-differential-equations for the correlation- and response-functions. 
Initial conditions are q{t,t) = 1 and r{t,t) = 1. It should be stressed that 
the above equations do not require equilibrium. They are therefore suited 
to deal with off equilibrium properties and aging. In general the correlation- 
and response-functions depend on both time arguments t and t' and not only 
on the difference t — t'. 

3.6 Equilibrium Dynamics in the Ergodic Phase 

Above the freezing temperature T^, in equilibrium, the correlation- and res- 
ponse-functions depend on the difference t—t' only and they obey fluctuation- 
dissipation-theorems, Eq.(jSZ|). The above equations of motion simplify con- 
siderably. The memory terms obey fluctuation-dissipation-theorems as well 

K{t) = -(3M{t) (52) 

and Eq. (j47j) reads 

(^^+7z)g(t) = /im + /?|M(t)-M(oo)g(oo)} - 13 j^ds M{t - s)q{s) (53) 

with M{t) = W'(^q(t)). Magnetisation m = m(t^oo) and fl = /i(t-^oo) are 
given by 

]Im = h + p {W'{1) -W'{q)}m (54) 
with q = g(t — » oo) , and 

JI= hm + T + (3 {W'{1) ~qW'{q)}. (55) 

For h = one has m = q = and JI = T + f3W'{l). The following 
discussion will be restricted to h = 0. The general case with h ^ is discussed 
in [221 and results will be shown later. 

For general t, Eq. (j53p can be solved numerically. Since we are dealing here 
with a purely dissipative system, q{t) < is required for all t. Replacing 
q{s) by q(t) in the integrand of Eq. lj^Hj) the following inequality is obtained 

< -m < q{t) - (3 W'{q{t)) (1 - g(t)) (56) 



or 



Tq{t)>W'{q{t)){l-q{t)). (57) 
Left and right hand side of Eq. (j37|) are shown in FigEl 



17 



A 1 q D 1 



Figure 9: Stability criterium Eq. ()57|l for h = 0. (A): p = 2 (B): p > 2. The 
left hand side is plotted for various temperatures. For T < Tc the stability 
criterium is violated. 

Assuming W{q) ~ for g ^ two cases have to be distinguished: 

A) p = 2: For h = a. second solution branches off for T < T^. The 
stability criterium is, however, violated for q{t) < q^. This means that 
a phase transition into an off equilibrium phase has to take place at 
T = Tf.. Since this solution branches off continuously from g = at 
Tc this transition is referred to as continuous transition. For h ^ the 
straight line representing the left hand side cuts the g = 0— axis at some 
negative value and there is no transition in finite field [HH] . 

B) p > 2: For T>Tc again a single solution q = exists. At Tc{h) a new 
solution qc fulfilling 

W'{qc) (l-gc)M = Te (58) 

shows up. This means that the resulting phase transition is discontin- 
uous. For T = Tc and q = qc the slope of the left and right hand side 
of Eq.()57|) has to be the same resulting in 

Tc = W"{qc){l-qc)-W'{qc). (59) 

The two Eqs. ()58|) and ()59|) determine Tc and qc- 




Figure 10: Phase diagram of the spherical p-spin interaction spin-glass [22] • 

This discussion can be extended to h ^ [22]. The resulting phase diagram 
is shown in FigHHIfor the p = 3-spin-glass with W{q) = ^q^. 
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For h < 0.41 the transition is discontinuous. The transition temperatures 
Tc{h) obtained from dynamics is higher than the one found in rephca theory. 
This surprising result might be explained by the proposal that the domi- 
nant states counted for in replica theory are not reached by the dynamics 
investigated (iOj . 

For h > 0.41 the transition is continuous and both theories give the same 
transition temperature. 

The results of a numerical integration of Eq. ()53|) for h = is shown in 
Fig^2 for temperatures slightly above and at the critical temperature T^. 
The asymptotic behaviour q{t) close to Tc can be analysed by identifying 
appropriate scaling functions for different regions in time and by matching 
nT\ 1^ . This resembles very much the crossover scaling analysis near 
multi critical points jUj. From Fig^J one can read off two characteristic 
temperature dependent time scales: A plateau time Tp(T) can be defined as 
q{Tp) = Qc and a second time scale Ta{T) characterising the final decay. It 
may be defined as q{Ta) = \qc- 



For t <C Tp{T) the correlation-function approaches some universal function 
q{t) qoit) as T ^ %. For K t < Tp{T) Eq.((SSl) is solved by 



This result is obtained by inserting Eq. ()60p into Eq. ()53|) with T = T^. For t— > 
oo the leading contribution is (9(1). The next to leading contributions 0{t~°') 
yield Eq. (j3Hj) and Eq. (j^ respectively. Collecting terms 0(t~^") Eq. (j^ is 
obtained. 

The scaling function ansatz for t ~ is 




1.0 



Figure 11: Correlation function q{t) for T > Tc |22j . 



qo{t) -^qc + Cot 
The dynamical critical exponent a is solution of 

r(l - a)' _ (1 - qc)W'"{qc) 
r(l-2a) 2W"{qc) 



(61) 



(60) 
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The contributions 0{T — Tc) to Eq. (p]|) for t Tp give 

Tpr~.(T-T,)-'/^\ (63) 
For Tp <^t <^ Ta Eq. ffSHj) is solved by 

qp{r)^~CpT'' (64) 

with b obeying 

T{l + bf {l-qc)W'"{q,) r(l-a)2 



r(l + 26) 2W"{qc) r(l-2a)' 



(65) 



This means that the two dynamical critical exponents a and b are not inde- 
pendent [inil221- 

The last scaling regime applies for t ^ Ta- The scaling ansatz is 

qit) = qait/ra). (66) 

Matching to Eq. (jU^ yields 

Qa{r) ^ qc~Cp for r (67) 

and 

Tp = r'J^'^^'^ (68) 

and with Eq.(|U^ 

r„ ~ (T-T,)-("+^)/2''^ (69) 

This means that only a single independent critical exponent exists. The 
scaling functions qo(t), qpir) and qair) have to be determined numerically. 
The imaginary part of the frequency dependent susceptibility is 

/'OO /'OO 

X"(cj) = / &tr{t) sm{ujt) = I3uj j dtq{t) cos{ujt) (70) 
Jo Jo 

where the second expression is derived from the FDT, Eq.(|37p. Near the 
main contributions are due to t ~ 1 and t ^ Ta, respectively. Inserting the 
scaling discussed above 

/•oo rcx> 

X"M = / dtqoit) cosiu;t)+P^I drg,(r) cos(^t) 

= X'>] + Xl{^j " (71) 

with 

a;, = l/r,~(T-T,)('^+^)/2'^^ (72) 
With Eq.dnni) for 

X'oioo) ^ l3coT{l - a) cos(|7r(l - a)) tu'' (73) 
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Figure 12: Susceptibility x"{^) computed from the numerical solution of 
Eq. (jKH|) foTh = and various temperatures characterised by TaiT), Eq. (jHn|) . 
For comparison the susceptibility for finite h in the region of continuous tran- 
sitions is also shown. 

is found. For u Ua the second contribution with Eq.(|Mj) gives 

x:{u;) ^ PcaT{l + b) cos(i7r(l + b)) u-\ (74) 
Combining both expressions x"!"^) shows a minimum at a frequency 

~ rl'^"^'^ ~ ~ (T - T,f'^\ (75) 

The complete susceptibility is plotted in FiglT^ The similarity with the 
experiments on glasses, shown in Fig|B^, is obvious, and indeed the data on 
CKN and other glasses have been fitted to the equivalent results of mode 
coupling theory [ini ESI 121] reproducing the data quite well. 

3.7 Off Equilibrium Dynamics and Aging in the QFDT- 
Phase 

The investigation of the dynamics for T < requires to specify some long 
time scale Too, and to consider the limit Too — ^ oo eventually. There have been 
several proposals for such a time scale: 

• For finite systems with N elements the equilibration time too is finite 
0. For large A^, Too ~ e^'^' g^. 

• Relaxing interactions J{t) J(t') ~ e^*~*')/^°° pUj . 

• Slow cooling with T(t) = (1 - t/xoo) |221 ES]. 

• Fast cooling, aging, with r^o = [HI US] . 

In the following we concentrate on aging. The system is rapidly cooled from 
high temperatures to some T < at t = 0. Alternatively a strong magnetic 
field could be switched off at t = 0. The requirement is, that the initial state 
is not correlated with the disorder. 
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Figure 13: Correlation function q{t + 1^, t^) and FDT-violation X{t + 1^, tw) 
for T < Tc and h = 0. The curves are plotted for fixed t + as functions of 
t/tiu. The range t/t ~ t ~1~ tiij corresponds to t^j ~ 1. 



Results from a numerical integration 133] of the dynamical mean field 
equations ()47ll5ip are shown in FiglT^ 

Because the system is not in equilibrium at T < Tc, correlation- and 
response-functions depend on both time variables t and t' , and not on t — t' 
only as in equilibrium. Furthermore fluctuation dissipation theorems, Eq.(|37|). 
are not fulfilled. A measure for FDT-violation is 

^M=T^^, (76) 

which is shown in the lower part of FigfTT^ This quantity can be relates to 
an effective temperature |46j 

T,s{t,t')=X~\t,t')T. (77) 

For short time X{t + t^,t^) ^ 1 is found. This means that the FDT is 
fulfilled, while q{t + t^,t^) decays from 1 to qc- A distance in phase space 
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can be defined as 



Dit,t^) = ^Jl(^{^i{t + t^)-^i{t^)Y^ = ^ l-q{t + t^,t^). (78) 

This distance grows until q(t + t^jt^) ~ Qc is reached, and the resulting 
Dc = a/1 — qc can be interpreted as the size of a typical valley in phase space. 
The system apparently equilibrates within such a valley, and escapes only 
after some characteristic time Tp(t^). This escape is irreversible, which is 
indicated by the deviation of X(t + t^,, t^) from 1 starting at t ~ Tp(t«,). 

For t > Tpitw) the parameter X{t + t^, t^) reaches a plateau at some Xc, 
which extends up to t ^ t^. In this regime a modified fluctuation dissipation 
theorem (QFDT) holds. A solution of this kind was first observed in the 
context of learning in neural networks exhibiting a discontinuous transition 
as well jni. It appears to be a common feature of discontinuous freezing 
transitions |ill22]. 

Investigating the asymptotic behaviour of the solutions of the dynamical 
mean field equations, Eqs.(jl3EII), the parameters q, JI, fn and Xc can be 
evaluated analytically P^. Especially qc is given by qc = PW'{qc) (1 — qc), 
which is of the form of Eq. ()58|) . indicating that the solution is marginally 
stable at all temperatures T < Tc. 

The crossover scaling analysis is similar to the one discussed for ergodic 
dynamics in the previous section, but the time scales Tp(tw) and Ta(tw) are 
now determined by the waiting timet^. 

Again dynamical critical exponents are introduced. They obey 

r(l-a)^ ^ ^ T{l + by ^ (1 - qc)W"'{qc) 
r(l-2a) 'r(l + 26) 2W"{qc) 

which is identical to the result Eq. ljU^ for T > Tc , except for the factor Xc. 
The behaviour of X{t + t^j, for 1 <C t -C and <C t <C Tq, respectively, 
is also ruled by power laws, as indicated in the lower part of FigHni The 
resulting exponents fulfil [121 US] 

a' = 3a + l 6' = 36-1. (80) 

Obviously for the existence of the QFDT-phase 6' > has to be fulfilled. 
b' = therefore marks a phase transition from the QFDT-phase (^-phase) 
to a phase with a hierarchy of long time scales (i3-phase) . Another criterium 
to be fulfilled is X(t,t') > 1. This gives rise to yet another phase (C-phase) 
which covers only a tiny portion of the phase diagram. 

The complete phase diagram is shown in FigfT^ Qualitatively the same 
phase diagram is found for the mean field model of a particle moving in a 
correlated random potential jl]. 

In phase A the plateau-time Tp and the characteristic time scale for aging 
Ta are again related by Eq. ljUHj) with the modified relationship, Eq. (f7^ . among 
the exponents. 
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Figure 14: Phase diagram of the spherical j9=3— spin interaction spin-glass 
|45j . O: Ergodic phase. A: QFDT-phase. i3 and C: Phases with a hierarchy 
of long time scales. The points a and b are discussed later. 



For t ty, it has been argued jlH [T2] that 

qit + t^,t^) =Q{h{t + t^)/h{t^)) (81) 

as long as X(t + tw, tw) ~ Xc. The function h(t) has to be monotonic, but is 
otherwise not determined. The arbitrariness in h{t) resembles the invariance 
postulated for the long time behaviour of the SK-Model jUJ [TU] . 

Q{x) can be computed numerically and Q{x) — qc ^ {x — I)'' for x 1. 
Especially for zero field the exponent 6=1 and 

Q{x) = ^ (82) 

X 

is found. The postulate of scaling with Ta{tw) = yields h{t) = t^. The 
numerical results [13] shown in Fig^J agree reasonably well with this form 
of scaling and in particular with c ~ 1/2. 

For the analogous case of the drift of a particle in a random potential j3] , 
Ta ~ , with /i = 1 — a(36 — 1)/ (2(a + 6) + a(36 — 1)), has been found in 
the ^-phase. In this problem the time scale = 1/ f drift plays the same role 
as for aging. Typical values are /i ~ 0.9 ■ ■ ■ 0.95. Scaling of the form given 
in Eq.(|66|) with Ta{t^) ~ is obtained from h{t) = e'^*''^^)'^/'' with r] = 1 — fx. 
The numerical solutions indicate r] > 0, but the longest time investigated 
is not long enough, to decide on the actual value of 77. 



3.8 Off Equilibrium Dynamics and Aging in the B- 
Phase 

We now turn to the i3-phase. FiglT^ shows correlation-function c(t + tw, t^) = 

q(t + t^,t^) — m(t + tw) m{tyj) and FDT-violation parameter X{t + t^j^tu,) 
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Figure 15: Correlation function c(t + 1^, t^) = q(t + 1^, t^) — m(t + 1^) 
and FDT-violation parameter X{t + t^, t^) for T = 0.5, h = 0.5 (^-phase) 
and h = 0.7 (i3-phase) 03j corresponding to point a and b in Fig|14l The 
curves are plotted for fixed t + . The range t/t ~ t -|- til, corresponds to 



for field and temperature values indicated as points a and b in FiglT^ The 
qualitative difference between A- and S-phase is clearly visible. 

In the ^-phase, for t ^ Tp and ^ 1, the correlation function approaches 
a scaling form 

C{t + t^, t^) Cait/Ta) with Ta ^ ty, , (83) 

whereas no such scaling is found in the ;B-phase. This is in some analogy to 
one step versus full replica symmetry breaking [S]. 

Within replica theory the p=3-spin-glass does, however, not exhibit a full 
replica symmetry broken phase at all. This is another difference in the phase 
diagrams obtained from dynamics and replica theory. For the same range of 
time, the FDT-violation parameter X{t + develops a plateau in the 

^-phase, consistent with the QFDT-solution, whereas no such plateau shows 
up in the i3-phase. 

3.9 Mean Field Dynamics of Systems with Short Ranged 
Interactions 

The slowing down of the dynamics near Tc in a disordered system is not 
expected to be associated with a diverging length scale, at least not on the 
level considered here. A diverging length scale is expected only for certain 
four point correlation functions, which are factorized in the dynamic mean 
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field theory or in mode coupling theory, e.g. 

C{Ri„t-t') = (y;,(t)y;,(t)y.,(f)</^,(t'))"'-(^.W<^.(i'))(</'iW</'i(^'))^- (84) 

Dynamic mean field equations for short ranged interactions can be derived 
adopting a factorization property 

%j{t, t') = fi,j ■ q{t, t') r,j(t, t') = fij ■ r{t, t'). (85) 

The corresponding factorization has been proposed and tested within mode 
coupling theory for supercooled liquids (THl HB] ■ The resulting equations for 
the time dependent parts, q{t,t') and r{t,t'), are of the same form as those 
derived for mean field models, Eqs. ()47ll^ . The definition of W{q), Eq. p^ . 
has to be modified, taking into account the short range nature of the inter- 
actions and the form factors fij. The crucial assumption is, to neglect four 
point correlations of the form given above. 

It should be mentioned, that a rather different picture, based on a droplet 
model description, has been developed [Hj . Both pictures, mean field models 
and the droplet model, reproduce certain features of short ranged spin-glasses, 
and no clear cut evidence in favour or against one of the theories seems to 
exist at present. 

3.10 j3-Spin-Glass with Relaxing Bonds and Cage Re- 
laxation in Supercooled Liquids 

In a supercooled liquid one may distinguish two kinds of motion. Each particle 
in a sense is surrounded by a cage formed by other particles. Within this cage 
it may vibrate with amplitudes much smaller than the interparticle spacing. 
This picture is of course simplified in the sense, that the cage is also vibrating, 
and each particle is also part of the cage of other particles. It is probably 
more appropriate to view this motion as a collective localized anharmonic 
vibration. The second type of motion is a rearrangement of the particles 
forming the cage, or an escape of a particle from its cage. This type of 
motion involves jumps over distances of the order of the interparticle spacing, 
and this is supposed to be an activated process [201 1211 • 

The mode coupling theory for supercooled glasses [121 [TB] deals with den- 
sity fiuctuations, and the resulting equations are identical to those derived for 
the p = 3-spin-glass in equilibrium, Eq. (j53p . The initial decay of the correla- 
tion function towards the plateau (see FiglTTj) is usually interpreted as being 
due to the motion within the cages, whereas the ultimate decay is supposed 
to result from the escape of a particle from its cage. On the other hand, acti- 
vated processes are not contained in the standard mode coupling theory. The 
extended theory [IH] captures those processes, and yields a rounding of the 
singularities found in the ideal mode coupling theory at Tc. This extended 
theory is still an equilibrium theory and can not describe aging [2^. 

In the following I sketch a slightly different picture |45, resulting in iden- 
tical mode couphng type of equations above T^.. This formulation allows. 
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however, to deal with off equihbrium properties. It is based on a distinction 
of the two types of motion described above. 

We neglect for a moment the activated process of rearrangement of the 
cages and consider the anharmonic vibrations only. The position of a particle 
at time t may be written as 

Xi{t) = Ri + Ui{t) (86) 

where 'Ui(t) is the momentary displacement of particle i from its mean po- 
sition Ri. The configuration of the Ri can be viewed as inherent state P^]. 
The dynamics of the displacements Ui{t) is treated within the selfconsistent 
anharmonic phonon theory ^H] , developed originally for quantum crystals or 
strongly anharmonic crystals at high temperatures. The main idea of this 
approach is, to replace the harmonic and anharmonic coupling constants by 
their averages over thermal or quantum fluctuations, e.g. 

{V)., = jdrV\r)g,,{r). (87) 

V{r) is the pair potential and gij{r) the static pair correlation function be- 
tween particles i and j. Employing the factorization property Eq. (jH3j) and 
retaining cubic anharmonicities, one ends up with equations identical to those 
of the dynamic mean field theory, Eqs. ()47ll^ . 

The coupling constants (^)j... are calculated for fixed mean positions Ri. 
These positions are, however, random and the coupling constants are therefore 
random variables as well. They play the role of the Jj... in Eq. ()18|l . Taking 
into account the activated motion of the Ri(t), the random couplings Ji...(t) 
now depend on time as well. Under the assumption that they are Gaussian 
distributed random variables, Eq. (fT!?|) holds in the modified form 

j~mln = ^{k, ■■■ + ••■ }w,Gjit - 1'). (88) 

Since the couplings J...{t) are dynamical variables, corresponding correlation 
functions have to be introduced 

JUt)J,...{t') = -r^Akj ■■■ + ■■■ W,Fj{t - t'). (89) 

The resulting mean field equations are unchanged, except for the modified 
memory terms, Eq. (j5Up . 

K{tX) = W"{q{tX))Gj{t-t')r{t,t') + W'{q{t,t'))Fj{t-t') 
M{t,t') = W'{q{t,t'))Gj{t-t'). (90) 

Assuming an activated process for the cage relaxation, 

Gj{t) = e-*/"=^^" and Fj{t) = -^q-'I^--^^ (91) 

''"cage 
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is chosen. 

The choice /?j = and Tcage ^ cc is a reahsation of quenched disorder 
[TUj . The results for 3> Tcage are similar to those obtained for aging in 
Sect ]3. 71 with Xcage replacing t^. This holds in particular for the algebraic 
decay of the correlation function q{t) — gc ~ 

The situation for (3j = (3 is quite different. This means that the dynamics 
of the cages (bonds) can also equilibrate, and that fluctuation-dissipation- 
theorems hold for ^ Tcage- Numerical integration of the dynamic mean 
field equations yields results shown in FiglT^ 
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Figure 16: Correlation-function q(t + t^, t^) for T < Tc, cage relaxation time 
T"cage = 10^ and various waiting times. For comparison the correlation function 
for quenched disorder (xcage = oo) is also shown (thin lines). 

There are now two long time scales, the waiting time and the cage relax- 
ation time Tcage- There is yet another characteristic time scale fcage = Tp (Tcage) 
associated with the cage relaxation time, where Tp is the plateau time intro- 
duced in Sect ]3. 71 

For t + t^i; < Tcage the results found for quenched disorder are recovered. 
This means that the approach and departure of the correlation-function from 
the plateau qc is ruled by power laws, fiuctuation-dissipation-theorems are 
violated for t > Tp{tyj) and aging is observed. 

For tyj ^ Tcage the system equilibrates, a new plateau value > qc is found. 
The ultimate decay of the correlation function is now ruled by Tcage- No 
critical behaviour, i.e. no power laws in t, are found. The exponential decay 
of the correlation function towards the new plateau value is associated with 
a time scale Tp(T) ~ (Tc — T)^^. 

The susceptibility x"{^)^ calculated for the p-spin-glass model with relax- 
ing bonds, depends on the assumptions made for Fj{t) or (3j. For equilibrated 
cage relaxation (3j = (3 a. "knee" shows up at Up ~ I/t^, and x"('^) ~ ^ for 
uj < ojp whereas x"('^) ~ ^'^ for ^ > ^p- Such a knee has actually been 
observed in CKN but this observation was later withdrawn. The knee is 
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Figure 17: Temperature dependent time scales for the p-spin-glass with re- 
laxing bonds. For the waiting time indicated in the figure the equilibrium 
and aging regimes are marked 

absent assuming non equilibrated cage dynamics with Fj = 0. 




Figure 18: Susceptibility for the p-spin-glass with relaxing bonds. 

"Quenched disorder" calculated with /3j = 0, "annealed disorder" calculated 
with /3j = /3. 

It has to be mentioned that the relaxation of the cage configurations 
has been put in by hand. In particular it is not infiuenced by the critical 
slowing down of the anharmonic vibrations. This is a consequence of the 
mean field model and the A^-dependent scaling of the interactions, Eq. (j20|l . 
More realistic models with short range interactions would yield a coupling 
of the dynamics of the interactions J and the anharmonic vibrations, even 
within the mean field approximation discussed in Sect J3. 91 This has, however, 
not been worked out. 

The picture outlined in this section attributes the dynamics of a super- 
cooled liquid, on a time scale shorter than the life-time of cages, to anharmonic 
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Figure 19: Coupled 1-2-phonon oscillation in an anharmonic potential. 

vibrations. The instability at is attributed to a softening of a coupled oscil- 
lation of the center of mass and the width of the thermal cloud, representing 
the motion of a particle in its cage. This motion is sketched in Fig^J The 
idea of a coupled 1-2-phonon process was originally developed for quantum 
crystals [IH]- The critical slowing down near Tc within this picture is not 
directly related to the glass transition, and the actual glass transition at Tg 
has to be attributed to a slowing down of the activated cage relaxation. 

The standard interpretation |2ll is, to identifying the dynamic 

transition temperature Tc with the crossover region in glasses, and the onset 
temperature of replica symmetry breaking Trsb with the Kauzmann temper- 
ature, which has to be lower than the glass transition temperature Tg. A 
comparison of T^ and Tg for various glass formers, as done in Fig|31 reveals 
that both temperatures can be rather far apart, whereas model calculations 
jnni result in differences between and T^isb niuch smaller than the values 
observed. 

4 Outlook 

In this lecture I have concentrated on the dynamics of disordered systems 
treated in mean field theory. I have discussed models within physics but also 
models for various questions outside physics. Among those are optimization 
problems, neural networks and problems related to economics. Especially 
for the applications outside physics, where spacial dimension is not relevant, 
mean field theory is in many cases exact. The dynamic mean field theory is 
quite rich because the order parameters are functions of two time variables, 
the correlation- and response-functions. At high temperature or strong noise 
level, the systems are ergodic. Lowering the temperature or noise level, a 
freezing transition shows up and the system becomes non ergodic. 

Even on the level of mean field dynamics there are open questions, for 
instance the arbitrariness in the time parametrization function h{t) discussed 
in Sect J3.7| or the scaling form of correlation- and response-functions in the 
aging regime of the i3-phase, discussed in Sect ]3. 81 A more thorough under- 
standing of the relationship between dynamics and replica theory is desired. 
In particular there is no obvious prescription of how to calculate entropy or 
free energy from dynamics. These quantities are, however, of prime interest 
in replica theory. 
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The coarse features of aging in spin-glasses are reasonably well described 
by dynamic mean field theory. There are, however, more elaborate experi- 
ments investigating special cooling and heating schedules ^Tj. It is not clear 
whether they can be understood on the basis of dynamic mean field theory. 

Glasses and spin-glasses are systems with short ranged interactions. Some 
of their properties are captured by dynamic mean field theory or replica the- 
ory. Spacial correlations are not taken into account and questions concerning 
the existence of diverging length scales can not be answered. On the other 
hand real space formulations, e.g. the droplet model for spin-glasses jl?] or 
models with kinetically constraint dynamics for glasses [^ZIIIHIj are concerned 
with spacial aspects and a better understanding of discrepancies and common 
features would be desired. 

After all, it is surprising that the dynamic mean field theory for disordered 
systems can be applied to so many systems in physics and in other disciplines. 
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